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ABSTRACT 

We used the Hubble Space Telescope WFCS near-infrared camera to image the host galaxies of a 
sample of eleven luminous, dust-reddened quasars at z ~ 2 - the peak epoch of black hole growth 
and star formation in the Universe - to test the merger-driven picture for the co-evolution of galaxies 
and their nuclear black holes. The red quasars come from the FIRST-I-2MASS red quasar survey and 
a newer, deeper, UKIDSS-I-FIRST sample. These dust-reddened quasars are the most intrinsically 
luminous quasars in the Universe at all redshifts, and may represent the dust-clearing transitional 
phase in the merger-driven black hole growth scenario. Probing the host galaxies in rest-frame 
visible light, the HST images reveal that 8/10 of these quasars have actively merging hosts, while 
one source is reddened by an intervening lower redshift galaxy along the line-of-sight. We study the 
morphological properties of the quasar hosts using parametric Sersic fits as well as the non-parametric 
estimators (Gini coefficient, M 20 and asymmetry). Their properties are heterogeneous but broadly 
consistent with the most extreme morphologies of local merging systems such as Ultraluminous 
Infrared galaxies. The red quasars have a luminosity range of log(Lboi) = 47.8 — 48.3 (erg s“^) 
and the merger fraction of their AGN hosts is consistent with merger-driven models of luminous 
AGN activity at z = 2, which supports the picture in which luminous quasars and galaxies co-evolve 
through major mergers that trigger both star formation and black hole growth. 

Subject headings: quasars 


1. INTRODUCTION 

Most galaxies in our Universe have supermassive black 
holes (SMBH) at their centers, which are thought to 
have grown in earlier epochs when galaxies had a higher 
gas content. Quasars — highly luminous evidence of 
rapidly accreting black holes — provide insight into 
this important stage of galaxy evolution. In particu¬ 
lar, dust-reddened quasars, such as those investigated in 
this project, are the most intrinsically luminous steadily- 
emitting objects in the Universe and may represent an 
intermediate stage between galaxy mergers and lumi¬ 
nous blue quasars, which eventually become quiescent 
SMBHs. 

Models of gas rich galaxy mergers have been proposed 
to explain the observed link between the growth of super- 
massive black holes (SMBHs) and their host galaxies by 
an evolutionary scenario in which the growing black hole 
moves from a heavily-enshrouded high-accretion phase 
as in observations of some local ultraluminous infrared 
galaxies (ULIRGs) through a brief “blow-out” phase 
where winds and outflows clear the obscuring dust, to 
an unobscured, blue quasar which later becomes a quies- 
cent black hole (Sanders et al.|1988 [Hopkins et al.|2d0^ 


20081. In this scenario, the so-called “blow-out” phase 
is expected to appear as a reddened Type I quasar (i.e., 
showing broad emission lines in its spectrum) with high 
Eddington ratios and strong outflows. Systems in this 
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phase are also expected to still show signs of the earlier, 
or possibly still ongoing, merger in the form of tidal tails 
and disturbed morphologies. 

Dust-reddened quasars (or “red quasars”), found pre¬ 
dominantly in samples of matched radio and near- 
infrared sources, satisfy the conditions for this blow-out 
phase. The largest sample of red quasars comes from 
the FIRST-2MASS sample (called F2M hereafter; Glik- 


il. 120041 [20071 

2012D . More recently, [Glikman et" al.j (|2013p matched 

the IIRST radio survey to the U nited Kingdom Infrared 
Telescope (UKIRT) Infrared Deep Sky S urvey (UKIDSS; 
Lawrence et al. 2007 Warren et al. 2007) over a pilot area 
of 190 deg^, identifying twelve new, fainter red quasars, 
three of them above z ~ 2 (prefixed with UKFS). In to¬ 
tal, our group has identified over 135 unique red quasars 
i n the redshift range 0. 1 < z < 3. 


Urrutia et al. 120091 IGlikman et al. 
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Glikman et al. (2012) showed that, when corrected for 
extinction, I'2iVI red quasars are more luminous than 
blue quasars at every redshift. In addition, the frac¬ 
tion of red quasars increases with intrinsic luminosity, 
in agreement with the blow-out model. Other surveys 
of obscured AGN (e.g., Assef et al. 2012) hnd that the 
obscured fraction rises with decreasing luminosity con¬ 
sistent with a higher covering fraction from a dusty, ax- 
isymmetr ic structure and consiste nt with a unified model 
of AGN (Urry & Padovani 1995). At higher luminosi¬ 
ties, this so-called “rece ding torus” is exp ected to vanish 
due to dust sublimation ( Lawrence|1991 ). Therefore, the 
reddening seen for red quasars is probably due not to a 
circumnuclear torus but to dust distributed throughout 
the host galaxy that is created during a merger-induced 
starburst. We are able to disentangle the two effects with 
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the F2M sample because F2M red quasars are selected 
to have broad emission lines in their spectra indicating a 
viewing angle that avoids toroidal extinction. 

A subsample of 13 F2M quasars at z ^ 0.7 were im¬ 
aged with the Hubble Space Telescope {HST) Adv anced 
Camera for Surveys (ACS) by|Urrutia et al.|((2008|. The 
images showed that F2M quasars are hosted by an un¬ 
usually high fraction of mergers and/or interacting sys¬ 
tems. While studies of the host galaxies of unreddened 


quasa.rs find a merger fraction of < 30% (Dunlop et al. 


2003 Floyd et al. 2004), the red quasar merger traction is 


S5%. in addition, S'pjteer observations of the same 13 red 
quasars reveal that most of these objects are accr eting at 
very high Eddington rates (Urrutia et al. 2012). These 
results support the picture that red quasars are an early 
dust-enshrouded phase of quasar-galaxy co-evolution. 

Despite the successes of these aforementioned merger- 
driven models, recent observations suggest a more com- 
plicated picture for AGN ge nera lly. Schawinski et al.| 
(2011), Simmons et al. (2012) and Kocevski et al. 


(I2012|) 

dected 


showed that most moderatAluminosity A-ray-sel 
AGN (10"^^ erg s“^ < Lx < 10'^^ erg s“^) at 1.5 < 2 < 3 
reside in undistur bed, disk-dominated galaxies. |Schaw-| 


inski et al. (2012) showed that this is also the case for 
Lboi ~ 10^^ erg s“^ heavily obscured quasars at z ^ 2. 

Although some theoretical models imply that the pres¬ 
ence of a disk does not in-itself eliminate the possibility 
of a merger, as disks can surviv e mergers or can be re¬ 
built quickly ( Puech et al.||2012 ), many models allow for 
stochastic accretion to dornmate at low luminosities (and 
low black hole masses and/or low accretion ra te) while 


mergers drive fueli ng at high luminosities (e.g., [Hopkins 


k. Hernquist 2006). This dependence on luminosity is 


supporte d by a recent meta-a nalysis of mergers in AGN 
hosts by Treister et al. (2012) in which the merger frac¬ 
tion rises monotonically over three orders of magnitude 
in bolometric luminosity, with the highest merger rate 
(85%) at the highe st luminosity bm (L ~ 10"^® erg s“^) 
represented by the Urrutia et al. (2008) HST imaging 
study of red quasars. 

However compelling, this picture has not been tested 
at z ~ 2, the epoch of peak quasar activity, where most 
SMBH growth is believed to occur. A handful of lu¬ 
minous blue quasars have been studied at 2 ~ 2 with 
the Near-Infrared Gamera and Multi-Object Spectrom¬ 


eter (NIGM OS) on HST (Kukula et al. 2001 


et al. 

2001 

mova 

of tl 


Ridgway 
jy the re- 


morphologies; only host luminosities were reported. This 
may be because obvious signs of mergers have already 
faded in these mature systems. In the present study, we 
use HST observations with the near-infrared detector of 
the Wide Field Gamera 3 (WFG3/IR) to examine red 
quasars in the rest-frame optical at z ~ 2, to probe this 
critical stage of black hole growth and galaxy evolution. 

The paper is organized as follows. Section describes 
the parent samples and selection of z ~ 2 red quasars ob¬ 
served with if5TWFG3/IR. Section [^describes the HST 
observations and the data reduction procedure. Section 
[4 provides photometric and source count statistics for the 
^served fields and Section [^describes multi-component 
parametric fitting to separate the quasars from their host 
galaxies. We discuss the derived properties of the quasars 


and their host galaxies in Section followed by a discus¬ 
sion of the individual quasars in Section We address 
the implications in Section and summarize our findings 
in Section [9l 

Throughout this work we quote magnitudes on the AB 
system, unless explicitly stated otherwise. When com¬ 
puting luminosities and any other cosmology-dependent 
quantities, we use the AGDM concordance cosmology: 
Hq — 70 km s“^ Mpc“^, VLm = 0.30, and Ha = 0.70. 


2. QUASAR SAMPLE SELECTION 

The F2M red quasars at z ~ 2 present an ideal sam¬ 
ple to test merger-driven co-evolution in the highest lu¬ 
minosity regime. Figure shows the de-reddened K- 
band absolute magnitude for red quasars (colored cir¬ 
cles) vs. redshift. For comparison, we plot with black 
dots unreddened blue quasars with FIRST det ections 
from the F IRST Bright Quasar Survey (FBQS; I Gregg 
eU^llQ^ and SDSS-LUKIDSS ( |Peth et al.|[20Tl| ). The 
quasars marked with thick circles are the t hirteen objects 
previo usly studied with ACS on HST by Urrutia et al. 
(2008). The sources marked with stars are the eleven 


quasars studied in this work. The sample spans about 
two magnitudes in intrinsic luminosity and the redshift 
range of 1.7 < z < 2.3. 

As Figure]^ shows, F2M quasars are among the most 
luminous quasars in the Universe after correcting for 
extinction. At the redshift range of our sample, the 
WFC3/IR 77—band is comparable to the ACS 7—band, 
which makes this study a high-redshift a n alog of the 
z ^ 0.7 sample studied by |Urrutia et al.| (|2008|). Al¬ 
though the z ~ 2 sample is more luminous by ~ 2 — 3 
magnitudes, this increase in luminosity is consistent with 
the ev olution of Ls, in th e quasar luminosity function 
(QLF; Groom et al. 2009) and thus means we are sam¬ 
pling tne same portion of the QLF in both studies. 

We selected for HST imag ing all ten F2M red quasars 
from Glikman et al. (2012) with 1.7 < z < 2.3 that 
had a near-infrared spectrum. To extend our sample to 
fainter magnitudes , we added an eleve nth quasar from 
the UKFS sample ( Glikman et al.|[2013 ) which also pos¬ 
sessed a near-infrared spectrum. Tablefdlists our targets, 
their near infrared magnitudes, their TF^and absorption 
and redshift. Note that these quasars experience 7F-band 
absorption ranging between Ak = 0.3 and Ak = 1.2 
magnitudes. 


3. OBSERVATIONS AND IMAGE REDUCTION 

The Hubble Space Telescope WFC3/IR camera is sen¬ 
sitive to near-infrared wavelengths from 800 to 1700 nm 
and offers 15 different filters of narrow, medium and 
wide wavelength transmission range. The pixel scale is 
0.12825"/pixel and the field of view of the IR detector 
is 136" X 123". For this study, we used the F160W, 
F125W, and F105W Liters, whose effective wavelengths 
are 1536.9 nm, 1248.6 nm, and 1055.2 nm, respectively. 
Images of each quasar were obtained in two filters - the 
F160W filter and either the F125W filter (for the three 
quasars with the highest redshifts) or the F105W filter. 

Figure shows the near-infrared spectra of the eight 
quasars imaged with F160W and F105W and the trans¬ 
mission curves of those filters. Figure shows the same 
but for the three highest redshift sources, whose blue im¬ 
ages were taken with the F125W filter. Note that in no 
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Figure 1. Dereddened K-band absolute magnitude of F2M and 
UKFS quasars vs. redshift. The respective iC-band limits of the 
F2M survey (if as < 17.4) and UKFS {Kj^g < 18.9) are indicated 
with dashed and long-dashed line. Colored circles correspond to the 
amount of reddening, as defined in the legend. The small black dots 
are blue, optic ally-selected quasa rs from the FIRST bright quasar 
survey (FBQS;|Gregg et al.|1996|l and radio-detected quasars from 
a deeper SDSoTUKlUbS catalog ( |Peth et al.||20lT| l for which no 
redden ing is assumed. The HST-imaged objects troin |Urrutia et al.| 
l|2008|l are marked with thick black circles and the sources studied 
in this work are marked with stars. 


case is the broad WFC3 F160W filter contaminated by 
the strong Hct emission line. The avoidance of strong 
emission lines and the large amounts of extinction in 
these quasars minimizes contamination from the quasar 
and allows for better sensitivity to low surface brightness 
features in the host galaxies. We also image our sources 
in a bluer filter to sample the host galaxy light below 
the 4000A break (marked with a vertical dashed line in 
the Figure), enabling us to study the host galaxies’ star 
formation properties. Table lists the WFC3/IR hlters 
used. 

The bottom panel of Figure shows the transmission 
curve of the three filters used in this study (solid lines) 
compared with the Y, J and Ks bands used in ground 


based imag ing, such as UKIDSS (dashed lines; Hewett 
et al.|[2006 ). 


We observed most of the quasars with both filters over 
a single orbit, reaching a 3a surface brightness of ^24 
AB magnitudes arcsec“^ per pixel-area in all bands. We 
report this value in a columns (11) and (14) of Table 
(see ij^for details on the surface brightness depth of the 
images). 

The WFC3/IR detector is capable of non-destructive 
readouts (NDRs) during an exposure and has several 
options for NDR sequencing depending on the dynamic 
range desired in an image. We observed our sources in 
MULTIACCUM mode using the STEP sampling which is 
a log-linear non-destructive readout mode that prevents 
saturation of bright stars and allows a broad dynamic 
range in a single exposure. Depending on the total expo¬ 
sure time, we used the STEP25, STEP50, or STEPlOO 
sequences. The observations were done in a 4-point box 
dither pattern, which helps improve the resolution of the 
final reduced images and enables sub-pixel sampling of 


the point spread function (PSF) which we drizzle to a 
final pixel scale Of 0.06"/pixel. The total exposure times 
of the reduced images are listed in columns (10) and (13) 
of Table [T for the blue and red filters, respectively. 

Eigures 4] and show the HST images of the eleven 
quasars. Visual inspection suggests that most have a 
nearby companion and/or a disrupted host. Although 
visual inspection has been used in pre vious works to 
ident i fy mergers in high redsh ift systems (Kocevski et al. 
2012| Schawinski et al. 2012) we also performed careful 
two-dimensional modeling of the point source plus host 
galaxy (presented in However, to better define the 
parameters of our images we first made astrometric and 
photometric measurements of sources in each field and 
used them to determine the significance of the nearby 
sources as true companions. 


4. ASTROMETRY, PHOTOMETRY AND SOURCE 
DISTRIBUTION 


Although the relative astrometry of the WEC3/IR im¬ 
ages is highly accurate (10 mas, according to the WEC3 
Data Handbook), the absolute astrometry of the images 
can be offset W as much as 1.5" from the absolute as¬ 
trometric gricj^ To determine the offsets and correct for 
them, we extracted source catalogs from each reduced 
science frame plus its associat ed weight image using the 
SExtractor software package (Bertin & Arnouts 1996) 
with a 5cr detection limit. We used these catalogs to 
match to SDSS for absolute astrometric correction and 
to UKIDSS to test and correct for deviations from the 
photometric zero-point provided by the WFC3 Guid^ 

We identified point sources by performing differen¬ 
tial photometry on our catalogued sources and plotting 
Amag in two apertures (in this case 12 and 20 pixels, or 
0.72"and 1.2") versus the MAG_AUTO which is a better 
estimate the source’s total flux than the aperture magni¬ 
tudes. An example of this analysis for the F160W filter is 
shown in Figure]^ Since SDSS matches will necessarily 
be brighter than most of the objects in the HST image, 
we restrict SExtractor to 5a sources brighter than 22.5 
magnitudes (AB) whose differential aperture photome¬ 
try lies along a constant locus, separate from galaxies. 
In each field we find between 6 and 20 stars that we an¬ 
alyze on a field-by-field basis to determine astrometric 
offsets. 

Using astrometrically-corrected images, we compared 
the position of the EIRST radi o image with the HST im¬ 
ages. Similar to the findings of Urrutia et al. (2008), we 
find that the radio peak overlaps the peak of the WEG3 
images. Since the angular resolution of EIRST images is 
5" (and pixel size of I'.'S), we could not compare struc¬ 
tural details of the radio source with any structure in our 
HST images. 

We also matched our catalogued sources to the 
UKIDSS DR9 LAS survey (using the multiple cone 
search tool in TOPGAT). Only five our fields have 
UKIDSS coverage (F2M134I, F2M1359, F2M1344, 
F2M0738, UKFS0030). The first four are imaged with 
FI05W while the only FI25W fields covered by UKIDSS 
is UKES0030. To perform a photometric comparison be¬ 
tween our WFC3 filters and traditional near-infrared fil- 


® http://www.stsci.edu/hst/wfcS/documents/handbooks/currentDHB/wfc3_( 
^ http: //www. stsci . edu/hst/wf c3/phot_zp_lbii 
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Figure 2. Near-infrared spectra of the eight targets imaged with F160W and F105W filters. The WFC3 filter curves are shown in each 
panel in blue and red, respectively. The locations of Ha and H/3-|-[0 III] are indicated by vertical dotted lines and the host galaxy’s 4000A 
break is shown with a dashed line. Note that in a few cases the Balmer lines are shifted into the atmospheric ab sorption bands and ar e 
not seen. For these objects we determined the redshifts from optical spectra, which were presented in Figure 6 of [Glikman et ^ 

The objects’ redshifts ensure that the strong emission from Ha does not enter the F160W bandpass, minimizing the quasar/galaxy ratio 
and enabling more accurate PSF subtraction. 


L F2M2222 + 0202 ? = 2.252 

F125w; 

F160W: 

4l<i 


i- F2M1531+2423 z = 2.287 

ilk 



j- UKFS0030 + 0025 ? = 2.299 

f 

Ik hvK 


: A' 

1 ' 1' 


F105W - - 

F125W - 

) \ F160W - 

' Y ■ ■ 

1 ' J 

H- 

1 


1.0x10* 1.5x10* 2.0x10* 

Observed Wovelength (A) 


Figure 3. Same as[^ except for the three highest redshift sources 
in our sample that were imaged with F160W and F125W filters. 
The bottom panel shows the transmission curves for the three 
WFC3/IR filters used in this study. For comparison we show the 
UKIDSS y, J, and H filters with dashed lines. 


ters such as those provided by UKIDSS, we shift the AB 
magnitudes of the point sources in our HST images to 
their Vega cou nterparts using the A B to Vega correc¬ 
tions given by |Hewett et al. (2006). We compare the 
F105W, F125W and h'16UW mE^mtudes to the UKIDSS 
Y, J, and H bands, respectively. We find that the two 
magnitudes are well fit by a line of slope unity. The in¬ 
tercept of the line fit represents an offset between the 
two magnitude systems (i.e., add 0.58, 0.49, 0.80 magni¬ 
tudes, respectively) which encompass differences in the 
filter transmission curves between the WFC3/IR broad 
band filters and the UKIDSS filters, as is evident in the 
bottom panel of Figure We report magnitudes in the 


rest of the paper using the WFC3/IR zeropoints on the 
AB system. The shifts are provided above to enable us 
to make quantitative statements in either system. 

We use the SExtractor catalogs to determine whether 
the excess of sources seen near our quasars is significant 
compared to the overall distribution of sources in each 
field. Figure shows histograms of the cumulative num¬ 
ber of matches as a function of source separation in C/.'S 
bins for the red quasars (red line) in the red and blue 
filters (excluding self-matches within 0"1). We compare 
this with the histogram of separations for the morpho¬ 
logically stellar sources (blue line) that we determined by 
differential photometry as described above. We exclude 
our quasars from the stellar source histogram. In addi¬ 
tion we measure the distribution of random coincidences 
(black line) shifting the positions of the stellar sources 
by 15" to the north and matching to the source cata¬ 
log. The mean size of the error for each set of matches 
is shown at the top right. 

We see that morphologically stellar sources show some 
excess matches compared to the random background, 
but that the red quasars show significant excess in both 
bands. A two-sided KS test comparing the distribution 
of sources near the quasars to the overall source distribu¬ 
tion (black line) as well as to the morphologically stellar 
sources in the F160W filter yields p-values of 2 x 10“® 
and 3 x 10“^° respectively, allowing us to reject the null 
hypothesis with greater than 5(T confidence. Further¬ 
more, KS test comparing the overall source distribution 
and the morphologically stellar distribution yields a p- 
value of 0.4, implying that for those distributions the 
null hypothesis cannot be rejected. Similar results are 
obtained for the blue filter data (I x 10“^, 3 x 10“® and 
0.6, respectively). The excess of sources within 4"5 is 
most pronounced. At z ^ 2 4''.'5 translates to a projected 
distance of ~ 38 kpc, and we adopt this distance as the 
upper limit for considering a companion system in our 
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Figure 4. Two color HST images of the 8 lower-redshift quasars studied in this paper imaged with F105W and F160W. Each row represents 
a separate object. The first column is the original image shown at a scale of S" x . The second column shows the residual image after 
subtracting only the point source component. The third column shows the model for all but the point-source component; the blank frame 
is a source to which no host component could be fit. The final panel shows the full residual including masked regions and is indicative of 
the overall goodness of fit. Evidence of me r gers a nd disrupted host galaxies is seen in most the sources. We apply the Red-Green-Blue 
color-combining algorithm of|Lupton et al.|(|2004|| to our images, and we average the count rate from the F105W and F160W images to 
produce the green frame. 
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Figure 5. Same asj^except here we show the three highest-redshift quasars imaged with F125W and F160W. 
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Figure 6. Differential aperture photometry for all the sources 
catalogued in our F160W fields showing the clear separation of 
point-like objects which lie in a tight locus (indicated by the blue 
region) from the cloud of extended sources. 

modeling. 


the peak positions of the chosen stars in the mosaic us¬ 
ing the daofind routine. To properly weight the stars 
in the PSF, we used the phot function to estimate their 
magnitudes using an aperture of 3.0 pixels and a sky an¬ 
nulus with an inner radius of 10.0 pixels and outer radius 
of 20.0 pixels. The pstselect algorithm then selected 
the brightest stars that were sufficiently separated from 
other bright stars, using a psfrad of 100 pixels and a 
fitrad of 4 pixels. The output from pstselect was 
then fed into the final PSF-making routine, psf, which 
verifies the chosen stars and fits a 0th order Gaussian 
to produce a hnal sampled PSF look-up table. Finally, 
seepsf task converts the lookup table to a FITS image 
of the PSF whose image size is 201 x 201 pixels. 

Figure shows the final PSF images for the F105W, 
F125W, and F160W filters, respectively. Tablej^lists the 
relevant parameters for each PSF, including the number 
of archival fields used to supplement our proprietary data 
in column (2), the number of stars used to create the PSF 
in column (3), the PSF’s full-width at half-maximum 
(FWHM) in column (4). 


5. IMAGE DECOMPOSITION 
5.1. Construction of the Point Spread Function 

In order to model the quasars, a point spread function 
(PSF) is needed to provide a standard of how real point 
sources (stars or quasar nuclei) are represented in a driz¬ 
zled image from a particular telescope. Previous studies 
of quasar host galaxies have devoted up to several orbits 
of HST time to obtain a deep, high signal-to-nois e ratio 


et al. 

2004 

Urrutia et al. 

2008 

). In more recent, sim- 

ilar w 

ork ( 

e-g-, 

Simmons et al. 

2012| Schawinski et al. 


2012 


2011), which examined the host galaxies of mod¬ 


erate luminosity X-ray-selected AGN in WFC3/IR data, 
the PSF was constructed from images of isolated stars 
across the field. These studies and our own examina¬ 
tion, found that the PSF does not vary strongly across 
the fielcj^ Therefore, we created one PSF for the each 
filter using stars in drizzled images of our own data as 
well as from archived observations in the same filter. We 
selected from the MAST archive all observations taken 
since 2012 January 1 with the WFC3 imager, in the IR 
aperture and the filters F105W, F125W, and F160W 
that used the same 4-point dither pattern. We chose 
only science observations and did not consider calibra¬ 
tion data. We visually inspected the results and selected 
for retrieval fields that showed a few isolated and bright, 
yet unsaturated, stars. We then processed the images 
in the same manner as our science data using the same 
astrodrizzle task, approximately doubling the size of 
our program data to create a high signal-to-noise PSF. 

We made 201 x 201 pixel cutouts around eligible stars 
and produced a single mosaic image of all the stars. We 
masked out any extraneous light sources and replaced the 
masked pixels with the approximated level of background 
noise so as to not interfere with the PSF fitting. 

We implemented the steps for PSF creation prescribed 
by the DAOPHOT package in the Image Reduction and 
Analysis Facility (IRAF) suite. We produced a list of 

® In addition, all our objects are observed within 10 pixels of 
each other in the same location on the chip. 


5.2. Multi-component fitting with GALFIT 

Our objective is to study the morphologies of red 
quasar host galaxies and to determine whether mergers 
play an important role in their triggering. To do this in 
a quantitative manner we carri ed out host/point-sourc e 
decomposition using GALFIT (Peng et al. 2002, 2010[ ), 
modeling our sources with a PSh' plus as many host- 
galaxy components as necessary to minimize the reduced 
while careful not to overfit the data. 

We model t he host gala xy components with a Sersic 
radial profile ( Sersic|[l9^ I given by the equation: 
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where re is the effective radius within which half the 
total flux is contained and the surface brightness at the 
effective radius is Eg. The parameter n is referred to 
as the ‘Sersic index’ and determines the concentration 
of the light profile. A prohle with n = 4 represents the 
light distribution of a classical bulge, while n = 1 is an 
exponential disk which hts a classical disk. When n = 0.5 
the function is a two-dimensional Gaussian profile. The 
parameter k is tied to n, so that the Sersic index is the 
sole parameter that independently determines the radial 
light profile. 

GALFIT fits a Sersic profile to an image by adjusting 
the following parameters and minimizing using the 
Levenberg-Marquardt algorithm: the x and y position of 
the profile’s center; the source’s integrated magnitude; 
Te', n; the axis ratio (b/a) and position angle. 

Previous quasar host studies hav e used a simily ap¬ 


proach with independent software ( McLure et al. 2000, 
Floyd et al.||2004 Urrutia et al.||20(J8| ) that hts separate 


PSb' and Sersic components in a two step process, where 
the PSF subtraction is performed initially by scaling the 
PSF to the peak flux in the quasar, followed by a host- 
galaxy fit. However, subtracting the point source hrst 
can bias the host galaxy fit, while GALFIT performs the 
point source plus additional component htting simulta¬ 
neously, reducing this bias. 
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Figure 7. Histogram of separations between red quasars and nearby sources detected in our SExtractor catalogs (red line). All three 
histograms are normalized to the number of input sources to represent a match probability. We compare these with the separation 
histogram for morphologically stellar sources in each field (blue line) and chance coincidence matches to a false catalog created by shifting 
the morphologically stellar sources to the north by 15”. We plot at the top left the mean size of an error bar that would be centered on 
a bin for each population. The distribution for sources in the F160W filter is shown on the left. Because of the smaller number of frames 
in each filter we combine the F105W sources with the F125W sources in the right hand panel. In both the red and the blue images, red 
quasars show a significant excess of companions within 4”5 of the quasar. 



Figure 8. Composite PSFs for the three filters used in our HST program. Left - F105W. Middle - F125W. Right - F160W. The images 
are shown with a logarithmic scaling with intensity ranges set to include 99.5% of the pixels and the image size is 201 X 201 pixels, 
corresponding to 12”6 X 12l'6. 


To begin fitting morphological components to our im¬ 
ages, we extracted a 201 x 201 pixel box centered on 
quasar’s peak, as shown in Figures]^ andSince GAL- 
FIT aims to fit all the flux in an image, it is important to 
mask out any additional sources of light or bad pixels in 
the image. We show the masked regions in the residual 
(4th) column of Figures]^ and 
Because reddening dims the quasar at shorter wave¬ 
lengths, the host-to-point-source flux ratio is larger in 
the F105W and F125W filters compared with the F160W 
filter. Therefore, to better decompose the two systems 
without the point source overwhelming the host flux, we 
performed the first fit for each source in the shorter wave¬ 
length filter. We then fit the F160W images indepen¬ 
dently, but informed by the results of the shorter wave¬ 
length fit^ 


Our ultimate goal is to remove the point source and 
recover flux from the underlying host galaxy, whose cen¬ 
tral position would be within 0"6 as the PSF, as well as 
any nearby companions. We define ‘nearby companion’ 
as any component with a separation between CK'6 and 
4'.'5 arcsec, which translates into a projected distance of 
^ 5 to ^ 40 kpc at the redshift of our quasars. How¬ 
ever, because an underlying host galaxy can be elusive, 
we began each fit with a single PSF plus Sersic index 
centered on any obvious companion galaxy component, 
plus a flat sky. If no obvious component was visible, we 
initially fit just a PSF plus sky. We then added an ad¬ 
ditional host galaxy Sersic component at the location of 
the quasar and compared the reduced value with the 
added component to decide whether an additional com¬ 
ponent was warranted. We adopted an added component 


® We experimented with fixing the parameters of the 
F105W/F125W components in the F160W images, but found that 


these often did not converge or yielded poorer results, with larger 
values. 
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if the reduced significantly improved (> 5cr in an 

F-test). 

In many cases, adding a Sersic component at the same 
location of the PSF in order to model the underlying 
host galaxy resulted in a significantly improved fit, but 
with GALFIT assigning unphysical parameters that im¬ 
ply a need to fit flux from an unresolved region at the 
position of the PSF (i.e., re < 3 pixels and n < 0.1). 
While unphysical, this added component improves the 

likely because it captures the residual noise from the 
single model PSF fit. However, since we are interested 
in decomposing the PSF from any discernible underlying 
host, in cases where an unphysical Sersic profile arose, we 
tried a fit with two PSF components whose position and 
magnitude GALFIT could adjust to account for residual 
flux from a single PSF fit. In all but two cases where a 
second PSF was attempted, the fit was significantly bet¬ 
ter (as determined by an F-test which gave cr > 10) and 
allowed for a more physical underlying host galaxy to be 
fit by an additional Sersic component. Although we use 
two PSF components to accommodate the flux from the 
quasar, we do not interpret this as evidence for a dual 
AGN, since the separations between the two components 
are smaller than the spatial resolution of the images and 
are more likely due to residuals. The difference in PSF 
magnitudes between the two components was generally 
smaller for the blue bands, ranging between 0.2 and 1.5 
magnitudes, compared with a range of 0.5 to 3 magni¬ 
tudes for the F160W images. 

To better understand the nature of the fits needing two 
PSF components, we tested this technique on ten bright 
stars found in five of the quasar fields. The brightness 
range of the stars spanned 6 magnitudes. The bright¬ 
est stars that we tested were two magnitudes brighter 
than the brightest quasar in our sample, and the faintest 
star that we tested was two magnitudes fainter than the 
faintest quasar (the quasars themselves span about 2 
magnitudes in brightness). We found, consistently, that 
single-PSF fits to the brighter stars resulted in large, 
symmetric residuals and high reduced values (> 20 
for the two brightest sources) which were improved sig¬ 
nificantly (by > 30%) with the addition of a second 
PSF component, but never succeeded at capturing all the 
flux. The fainter stars generally yielded good fits (typ¬ 
ical x^ < 3) with a single PSF and were not improved 
significantly with the addition of a second component. 

We combine the fluxes from the two PSF components 
into a single PSF magnitude in the following way: 

ruTot = -2.51ogio(10”°'‘'"^'’"’= -k 10-O-4 '"psf2)_ (2) 

In the two cases where an unphysical Sersic profile best 
fits the image, we ascribe the flux in the Sersic compo¬ 
nent to be part of the point source, combining the PSF 
magnitude (mpsp) and Sersic magnitude (ms) provided 
by GALFIT similarly: 

rriTot = -2.51ogio(10"°-^'"^s"' -k (3) 

We report these combined magnitudes as the PSF magni¬ 
tude in Table and indicate the origin of the magnitude 
with a footnote in the table. We do not show the second 
components in the galaxy model represented in the third 
column of Figures and 

In many cases, additional Sersic components are 


needed to better fit the companion systems, and those 
components may not have the same parameters, or even 
location in the two bands. We interpret these differences 
as potential regions of star formation which are bright 
below 4000A ( in the rest frame) showing up in the blue 
band, but not in the red band, or luminous regions of 
dusty star formation that would appear red rather than 
blue. While we can extract physical information about 
the companion systems from their Galfit parameters, in 
many of the sources the Sersic index that provides the 
best fit may not be physically meaningful (i.e., n < 1 or 
n > 4). Likewise the effective radius can range as high 
as 200". 

From Figures|^and[^we see that all but two (F2M0943 
and F2M2222) sources show nearby companions to the 
central quasar with disturbed mo rphologies. Another 
source (F2M1359, see Section 7.21 is fit by a relatively 
smooth and centrally located galaxy and may also not 
be hosted by a merger or even intrinsically reddened. 

5.3. Uncertainties in the multi-component fits 

Although the Sersic profile may provide some physi¬ 
cal insight into the light distribution of a galaxy, there 
can be multiple components and tidal tails that are not 
well approximated by the shape of the profile leading to 
large errors in the fitting parameters. Here we discuss the 
uncertainties of our results, keeping in mind that the un¬ 
certainties reported by GALFIT represent only the errors 
from the covariance matrix and do not generally account 
for other sources of error; reporting only these sources 
of uncertainty may underestimate the true errors. How¬ 
ever, there have been multiple studies characterizing the 
additional uncertainties in multi-component galaxy fits 
in the presence of a bright AG N. The first and largest of 
these studies was performed by Simmons & Urry (20081 
who simulated more than 50,0(30 images of AGiN and host 
g alaxies with a range of gal axy and nuclear properties. 


Simmons & Urry (20081 conclude that recovery of 
AGJN and host galaxy characteristics is generally very re¬ 
liable. AGN and galaxy properties are accurately recov¬ 
ered even in automated batch-mode fits, including cases 
where the host galaxy is outshined by the nuclear point 
source. In individually fitted images where the specifics 
of each source and image may be properly accounted 
for (e.g., companion galaxies and stars, noise variances 
across the image), the uncertainties are reduced from the 
batch-mode case. 

The uncertainties in component Sersic indices are also 
reduced when the centroids of the bright point source 
are separated from those of the galaxy components by 
at least the FWHM of the PSF, as is the case for all 
the comp anion systems in this wor k. Because the error 


tables in Simmons & Urry 
coincident components and 


(20081 assume positionally 
jatch-mode fitting, their un¬ 


certainty values are high compared to what we expect 
for the true uncertainties in this work. Nevertheless, the 
predicted conservative uncertainties are useful as guide¬ 
lines here. The error tables predict typical uncertainties 
in AGN brightness of 0.25 mag, typical errors in host 
brightness of 0.5 mag, and typical errors in host Sersic 
index of (5„ ~ 0.7. In most cases this does not affect the 
assessment of a galaxy component as likely to be disk-like 
or bulge-like. Host galaxy sizes (re) are more uncertain 
and the uncertainty depends somewhat on the fitted mor- 












10 


Glikman et al. 


phology. We list in Table [^uncertainties for the underly¬ 
ing host parameters as a combination of reported errors 
by GALFIT and estimated additional uncertainties due 
to multi-comp onent AGN and host ga laxy hts from the 
simulations of Simmons & Urry (2008) added in quadra¬ 
ture. The conipanion systems are tar enough away from 
the PSF that we do not expect the AGN to affect their 
fitting errors and we list for them the uncertainties that 
are reported by GAL FIT. 


The error tables in Simmons & Urry (2008) assume 
no detailed follow-up from batch-mode nttmg is per¬ 
formed and estimate maximum errors due to position- 
ally coincident AGN and host galaxy centroids, so their 
use in this study produces errors estimates that are con¬ 
servative with respect to the true uncer t aintie s. While 
an analogous study to Simmons & Urry (2008) basedon 
WFG3/IR data at our depth and AGIN to host ratios 
would be the ideal method for estimating our uncertain¬ 
ties, such an effort is significant and beyond the scope of 
this paper. We use the Sersic components primarily as a 
means of capturing all the flux in the images so that we 
can study the colors and luminosities of the merging com¬ 
ponents to better understand the co-evolution of m ergin g 
galaxies hosting luminous red quasars (see Section 6.2). 


6. RESULTS 

With the fully reduced, PSF-subtracted and residual 
images in hand, we can study their surface brightness 
distribution as a function of distance from the central 
point source. Bec ause t he red quasars’ images are highly 
asymmetric (see (|6.2.1 ) the single radial profile that is 
typically plotted tor galaxy light distributions is insuf¬ 
ficient to describe the profile of these red quasar hosts. 
We plot in Figuresj^and [l0| the surface brightness distri¬ 
bution, /i, in mag arcsec^'^for all the pixels in all three 
images as a function of radial distance from the quasar’s 
peak emission in each filter. The gray shaded regions 
and black contours represent the flux in the original im¬ 
age, while the blue contours are the fluxes from pixels in 
the PSF-subtracted imagj^ We plot the residual im¬ 
age’s flux in orange, and use their values to determine 
the statistical limit of our observations. We compute 
the standard deviation of the flux for all pixels in the 
masked residual image (as shown in the fourth column 
of Figure [^ using the IDL procedure mimn which is part 
of the astronomy routine library. The horizontal dashed 
line represents this Scr limit and is quoted in Table [l] For 
comparison, we plot in Figure im the same for the PSF 
profiles. 

In all the panels, except for F2M0943, we see significant 
structure beyond the point source. F2M0943 is the one 
object that does not favor a merger. Its profile in both 
filters are similar to the PSF profiles, as shown in Figure 

m 

The dynamic range of our images, which we define as 
the ratio of the peak flux in the point source to a 5- 
tr detection threshold {D = Npsp/N^a), is between 960 
and 4000 in the red (FI60W) band. In the blue bands 
(FI05W and FI25W) the dynamic range is between 360 
and 4500. The large dynamic range of our images poses 

By ‘PSF-subtracted’ we mean the removal of all centrally 
concentrated light, including from a second PSF component or a 
concentrated Sersic component (e.g., F2M1341) 


the biggest challenge for our ability to study the low sur¬ 
face brightness host galaxy light after point source sub¬ 
traction. This is because the poisson noise in the point 
source goes as V-^psf, which implies that the residu¬ 
als after PSF fitting and subtract ing is on the order of 
y/Npsp. When D is large, then V-^psf >> which 
ranges between 160 and 370 for the red filter and 60 to 
340 in the blue filter. Residuals tens to hundreds of times 
brighter than the faint host galaxy features means that 
even with perfect PSF modeling the Poisson errors will 
dominate over the underlying galaxy. In our sample, the 
source with the largest dynamic range in both filters is 
F2M0943. 

Most of the systems studied in this sample show evi¬ 
dence for mergers hosting the red quasars, although their 
details are heterogeneous and complex. We note that in 
the redshift range of our sample (z = 1.7 — 2.3) man y of 
the tidal features seen in the z ~ 0.7 sample from |Ur-| 
rutia et al. ( 2008) wou ld disappear. This was showii by 


Schawinski et al. (2012 ) who simulated the appearance of 
four F2iVI quasar iroml Urrutia et al. (2008) in the WFC3 
H band when redshifted to z = 2. Nevertheless, some 
merger signatures are still evident in the images. We ap¬ 
proach the interpretation of our images below with this 
result in mind, cautious not to over-interpret the details 
of the fitted parameters. 

6.1. Red Quasar Properties 


Urrutia et al. (2008) found that once the point source 
was separated from tne host galaxy, the quasars them¬ 
selves (i.e., the fitted PSFs in that sample) had redder 
colors than the low-resolution photometry reported for 
the systems as a whole. In fact, the shift in color for 
these objects was larger for sources with redder total 
E{B — V) values (as determined from spectral fitting). 
They attribute this effect to an excess of blue light in 
the F475W filter, which is bluer than rest-frame 4000A 


break. Urrutia et al. ([2008 ) interpret this as an excess of 
young stars adding blue light from the host to the total 
integrated colors. 

Do we see a similar situation at 2 ; ^ 2, where the 
quasars are more luminous by ~ 2 — 3 magnitudes, on av¬ 
erage, than the 2 : ~ 0.7 sample? We use the SExtractor 
catalogs described in Section as well as the magnitudes 
from the Galfit modeling to examine the effect of sepa¬ 
rating the host galaxy light on the PSF colors. In six 
cases we detect the underlying quasar host galaxy and 
detect offset components in the other three cases. Our 
experimental design was int ended to measure the same 
rest-frame emission as in the Urrutia et al. (2008) work, 
with filters chosen to straddle the 4000A break. The 
FI60W filter corresponds to ^ 4850 — 5890A in the rest 
frame, depending on the redshift, while the FI05W fil¬ 
ter corresponds to ~ 3250 — 3870A and the FI25W filter 
corresponds to ~ 3790 — 3840A. We therefore discuss 
the colors of the quasars in terms of rest-frame U — V 
corresponding to the FI60W filter and either FI05W or 
FI25W. 

We plot in Figure the rest-frame U — V color of 
our quasars as a function of their reddening, E{B — V) 
from|Glikman et al. (2012), comparing the total color of 


the source, as measured by the the MAG_AUTD parameter 
in SExtractor, to the magnitudes returned by the Galfit 
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Figure 9. Shaded contours show the surface brightness in each pixel as a function of distance from the point source (radius) in arc seconds 
for the eight quasars in Figure The panel on the left for each object plots F160W data while F105W is shown on the right. Grey 
shaded regions and black contours show the flux from the unaltered quasar image, while the blue regions show the flux distribution from 
the PSF-subtracted image, and the orange regions represent the data from the masked residual images. The 3cr threshold is shown with a 
dashed horizontal line. In each panel, with the exception of F2M0943, significant structure (i.e., companion galaxies or tidal features) is 
seen immediately beyond the PSF. 



Figure 10. Surface brightness as a function of radius in arc sec¬ 
onds the three quasars shown in Figure!^ The panel on the left 
for each object plots F160W data while ul25W is shown on the 
right. Contours are colored the same as in Figure]^ 


modeling. SExtractor’s MAG_AUT0 parameter is the equiv¬ 
alent of a “total magnitude” encompassing > 90% of the 
total contiguous flux centered on the peak of the light 
profile and therefore includes the quasar plus the host. 
We plot the U — V color from the SExtractor photometry 
with black circles. Red diamonds show the PSF compo¬ 
nents’ colors, while green triangles show the combined 
PSF component plus Sersic component for sources that 
are better-fit by such an added component. The dotted 
lines connect the values for a given source to guide the 
eye. 

In general, quasars with higher E{B — V) values also 
have larger U — V colors. Although in most cases (6/11) 
the PSF-only component does have redder colors than 
the total MAG_AUT0 magnitude or the PSF plu s Sersic 
component, we do not see the clear trends that |Urrutia| 
et al. (20081 noted in the z ~ 0.7 sample. The primary 


PSl' component is redder than the combined PSFs or 



Rodius (orcsec) 

Figure 11. Surface brightness as a function of radius in arc sec¬ 
onds for the three PSF images (Figure]^. 


PSF plus Sersic model and may contain some blue host 
flux as was found in the z ^ 0.7 sample. However, be¬ 
cause the dynamic range is higher in the z ~ 2 sam¬ 
ple, compared with the z ^ 0.7 sample, and since the 
WFC3/IR spatial resolution is lower than the ACS reso¬ 
lution, separating the point source from the host galaxy 
for these systems is more challenging; we cannot say with 
certainty whether or how much of the additional blue flux 
in the added components is due to the quasar or young 
stars in the host. 

Consistent with the notion that merger-induced star 
formation adds blue light to the host, we note that the 
two sources with the smallest color difference between 
the different photometric measurements, F2M0943 and 
F2M2222, do not appear to have nearby companions or 
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Figure 12. We plot rest-frame U — V color of red quasars versus 
E{B — V) to investigate whether the removal of host galaxy light 
results in a redder quasar com ponent, as was found in the 2 ~ 0.7 
sample of|Urrutia et al.|(|2008|l. The filled circles are from the total 
magnitudes estimated by biiixtractor’s MAG_AUT0 (offset to the right 
by 0.005 magnitudes, for clarity), which we compare with the PSF 
magnitudes determined by Galfit (red diamonds). In sources that 
are better-fit by a second PSF component we plot with orange 
triangles the combined flux as computed by Equation While 
sources needing an added Sersic component to absorb residual PSF 
flux are shown with green triangles computed by Equation]^ When 
taken alone, the PSF component is typically redder than the to¬ 
tal source magnitude, as well as the combined magnitude from an 
added component consistent with the lower redshift sample. This 
may indicate the presence of blue light added by a nuclear star- 
burst, but could also be an indication of the inadequacies of the 
single PSF fits. Two quasars (F2M1531 and UKFS0030) have the 
same reddening, with E{B — V) = 0.32. To distinguish between 
them, we plot UKFS0030 with enlarged symbols. 


be actively merging. Table lists the magnitudes and 
colors of t he quasars using the three metrics shown in 
Figure 12 (see (Efor details on individual sources). 


6.2. Host Galaxy Properties 

Having separated the quasar emission from the under¬ 
lying galaxy light, we can explore some of the host galaxy 
properties and compare them to what is known about 
normal and star-forming galaxies at similar redshifts. We 
add up all the flux that Galfit assigns to the non-quasar 
Sersic components by summing the magnitudes in flux 
space, in a similar manner used for the point sources. 
We compute the total magnitude of the host plus com¬ 
panion galaxy in each filter following, 

Ns^r 

ruTot =-2.51ogio(^ (4) 


where Nser represents the number of Sersic components 
in our best fit. 

Our chosen WFC3/IR filters closely correspond to rest- 
frame U and V bands at z = 2. This is illustrated in 
Figure [T^ where we plot three rest- f rame galaxy spec¬ 
tral templates from Kinney et al.| (19961: Elliptical, 
Sc and a moderately reddened starburst template with 
0.39 < E{B — V)< 0.50, chosen to span a wide range in 
star-formation rates. We plot the Johnson U and V fil¬ 
ter curves in dashed purple and green lines, respectively. 


Figure 13. Comparison of three galaxy spectra from |Kinney| 
|et al.|(|1996f, namely, elliptical (red), Sc spiral (black), and star- 
burst (blue) to the WFC3 F105W (solid purple line) and F160W 
(solid green line) filters at 2 = 2, which correspond to rest-frame 
U (dashed purple line) and V (dashed green line) filters. 


showing the sampling of the SEDs at these wavelengths. 
We then shift the F105W and F160W bandpasses to the 
rest-frame at z = 2 and plot them with solid purple 
and green lines, respectively. The same can be done for 
the three highest redshift sources, where the F125W fil¬ 
ter corresponds to the U band. It is evident from this 
figure that the observed IR and rest-frame UV/optical 
curves overlap significantly, allowing us to compare the 
red quasar host colors with the U — V colors of galaxies 
across the Hubble sequence at comparable redshifts from 
the literature. 

In Eigure we plot the observed infrared colors of 
the red quasars as a proxy for rest-frame U — V versus 
F160W magnitudes (corresponding to rest-frame abso¬ 
lute V-band magnitude) of the red quasar hosts (green 
circles) and companions (orange circles). Eor a com¬ 
parison sample, we turn to the Cosmic Assembly Near¬ 


infrared Extragalactic Legacy Survey (C AND ELS; Gro- 


gin et al. 2011 Koekemoer et al. 2011) which consists 
of deep, multi-cycle observations of well-studied legacy 
fields with the HST WFC3/IR camera using the same 


three filters as in this work. Bell et al. (2012) and 


Lee et al. (2013) have studied the morphological and 


star forming properties of CANDELS galaxies out to 
z > 2 and provide a useful comparison set to our red 
quasar host galaxies. We plot with plus signs CAN- 
DELS galaxies that had matches to sources in the pub¬ 
licly release d catalogs of the UKID SS Ultra-deep Survey 
(UDS; from |Galametz et al.||2013|) and whose photomet¬ 
ric redshifts are between 1.7 < -^ < 2.3 an d stellar mass 
M* > 3x IO^^Mq, which Bell et al. (2012) states as their 
completeness limit. In addition, we plot with blue aster¬ 
isks CANDELS photometry in the GOODS-South field 
of AG N in the same redshift range from |Simmons et al.| 
(2012). We see that the AGN and normal galaxies lie in 
the same part of this parameter space. 

The nearby companions to the red quasars (i.e., the 
orange circles) with the smallest photometric uncertain¬ 
ties lie in the locus of CANDELS galaxy colors. How- 
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Figure 14. Left - We plot the observed F105W—F160W (or F125W—F160W when applicable) as a proxy for rest-frame U — V color 
versus the observed F160W magnitude of red quasar host galaxies (green circles) and companio n systems forange circles). For comparison 
we plot the same for galaxies between 1.7 < z < 2.3 from the CANDELS survey (black crosses |Bell et al.||2012|l, while blue asterisks are 
AGN in the GOODS field from|Simmons et al.|(|2012|l. Right - F105W—F160W (or F125W— FlbUW) versus rest-frame F160W magnitude 
determined by summing all the components associated with the quasar (green circles). Black crosses and blue asterisks have the same 
meaning as in the left hand panel. Open triangles are the total magnitudes of the sources (MAG_AUT0), plotted for comparison, to indicate 
the extent to which imperfect point source modeling may affect the galaxy colors. 


ever, four out of the five host systems (i.e., the green 
circles) in the left hand panel are among the most lu¬ 
minous galaxies. The magnitudes of these components 
may be contaminated by imperfect point source modeling 
and therefore their luminosities may be overestimated. 
However, since their magnitudes are computed from the 
model fits, our estimates in most cases exclude residual 
host galaxy light in, e.g., clumps of star formation, which 
the smooth Sersic fit does not account for, which could 
lead to a lowering of their estimated luminosities. 

The very high luminosities of the red quasar hosts is 
not unexpected given the relative volumes probed by the 
CANDELS and F2M/UKFS surveys. CANDELS is a 
very small-volume survey covering a total area of ~ 800 
arcmin^, while FIRST, 2MASS and UKIDSS are all very 
large-volume surveys, enabling us to find the most lu¬ 
minous sources such as these red quasars. Therefore 
we expect our objects to lie in the high-mass/luminosity 
side of the diagram. Since stellar masses have been com¬ 
puted for the CANDELS sources, we examine their stel¬ 
lar masses in the magnitude range in which our systems 
lie and find that they are all comparatively high mass 
systems, with M* > 3 x IO^^Mq. 

In the right hand panel, we plot the total magnitude of 
all the light from host plus companions for a given quasar 
(computed using equation]^ to examine what the total 
magnitude of the system might be under the assumption 
that all the components are part of a merger that will 
eventually coalesce. Here, the colors of five of our red 
quasars appear to continue along the locus of CANDELS 
galaxy colors, even when including light from the nearby 
companions. 

The colors of the F2M quasars span a very broad range, 
broader than the CANDELS galaxies, which is indica¬ 
tive of the complex nature of these putatively merging 
systems, where dust extinction (leading to red colors) 


competes with star formation (leading to blue colors). 
If the CANDELS galaxies are undergoing secular evo¬ 
lution, their colors ought to change more gradually with 
mass. We know there is significant dust in the red quasar 
systems, since their intrinsically very luminous quasars 
are dust-reddened. For some of the companions with 
very blue colors there may be unobscured star forma¬ 
tion, making the hosts extremely blue. In other hosts, 
the star formation is behind dust, making them very red. 

To check whether imperfect PSF modeling affects the 
colors of the host galaxies, we plot with triangles the 
uncorrected total magnitudes (MAG_AUT0 from the SEx- 
tractor catalogs) of the red quasars in the left hand panel 
of Figure The total colors of our sources are bluer 
than the colors of the galaxies, implying that the blue 
colors of their host galaxies might be contaminated by 
some contribution from the point source. 


6.2.1. Non-parametric Galaxy Properties 

Our parametric approach of htting Sersic profiles was 
largely intended for PSF subtraction and capturing the 
residual flux for separate analysis. We are cautious not 
to over-interpret the fitted parameters themselves since 
parametric fitting methods are often insufficient to de¬ 
scribe merging and irregular galaxies and fail to account 
for all the structure that is seen. Non-parametric htting 
techniques have therefore been develop ed to asses the 
degree of irregularity i n galaxy images (Abraham et al. 
2003 Lotz et al. 20041. The Gini coefficient, G - orig¬ 
inally an econometric tool devised to asses the wealth 
inequality of a population - has been shown to corre¬ 
late with other morphological parameters, such as the 
concentration index. Likewise, M 20 , the second order 
moment of the 20% brightest pixels in the galaxy, and 
the asymmetry, A, which depends on the residuals of an 
image after subtracting a 180 degree rotation from itself 
are used to morphologically classify galaxies and identify 
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mergers. 

In the following analysis we apply these nonparamet- 
ric measures to the PSF-subtracted images (second col¬ 
umn of Figure 13 for the blue and red filters sepa¬ 
rately. To mitigate the effects of residual flux from PSF- 
subtraction, we exclude the flux from a 4-pixel-radius 
circle around the quasar’s position. 

We chose not to consider two other nonparametric 
quantities often used to asses the degree of merging in 
galaxies: the concentration parameter, C, assesses the 
flux distribution in a source by the ratio of flux in a circle 
(ellipse) with a small radius (semimajor axis) to a larger 
one encompassing most of the flux; and the smooth¬ 
ness parameter, 5, measures dumpiness in a galaxy by 
subtracting a smoothed image from the raw to quantify 
small-scale fluctuations ( Conselice|2003 together with A 
comprise the “CAS” system;), i'fie concentration index 
depends on inner galaxy light which we cannot accurately 
assess; the Gini coefficient does a better job of measur¬ 
ing concentration for our sources. The smoothness pa¬ 
rameter loses efficacy at high redshift where a resolu¬ 
tion element (PSFfwhm 0.15 — 0.2") corresponds to 
^1.5 kpc. Therefore, in this paper, we concentrate on 
G, M 20 , and A to study the morphologies of red quasar 
host galaxies. 

Since the depths of our 22 images vary and since the 
redshift of our sample introduces surface brightness dim- 
ming as (1 + z)^ , we follow the technique outlined in 


Lotz et al. (20041 to generate segmentation maps that 
dehne a set of pixels to analyze that are above a uniform 
threshold for the entire sample. We compute the surface 
brightness per pixel, as a function of radius (where the 
center of the map is the peak position of the quasar), 
/r(r). We compute the Petrosian radius, Vp, defined as 
the radius at which the surface brightness is 20% the 
mean surface brightness within that radius, i.e.. 


?7 = 


Krp) 

fl{r < rp) 


with r] = 0.2. 


(5) 


The segmentation map contains all the pixels with /r > 
fl{rp). 

We then create an array, Xi , sorted in order of increas¬ 
ing pixel value, and compute the Gini coefficient using 
the algorithm: 


G = 


1 


Xn{n — 1) 


^(2^-r^-l)X„ 


( 6 ) 


where n is the total number of pixels in the segmentation 
map and X is the mea n of all Xi values. 


Urrutia et al. (20081 found a correlation between the 
Gini coefficient and hl{B — V), indicating that redder 
sources are more disturbed. The same plot for our 2 ~ 2 
red quasars (Figure left) shows large scatter and no 
such correlation. 

The Gini coefficient alone does not determine the de¬ 
gree of merging/disruption, since a de Vaucouleurs profile 
(n = 4) is more centrally concentr ated than a disk , yet 
is still a smooth light distribution. Lee et al. (2013) find 
that passive, elliptical, GANDELS galaxies at z ~ 2 have 
G between 0.4 and 0.7, with a higher mean than the star 
forming systems (0.53 versus 0.43) largely due to their 
light profile being concentrated in a single central peak. 


However, a merger with bright clumps of star formation 
will also have large G values, making this parameter by 
itself insufficient for identifying merging systems. 

The distinction between a centrally concentrated light 
profile and a clumpy merger can be made when G is 
combined with a second parameter, such as M 20 , which 
represents the second order mom ent of the 20% bri ghtest 
pixels in a galaxy’s light profile (Lotz et al. 2004). M 20 
is defined as follows: 


M 20 = log 10 


Altat 


where 


Mu 




(7) 


( 8 ) 


Here n and Xi are the total number of pixels in the seg¬ 
mentation map and the flux per pixel, as defined for 
equation [3 We sum over Mi while Yi^i < 0.2Atot- 
Lotz et alT ( 2004 ) showed that smooth light profiles cor- 
respond to low values of M 20 , while extended sources 
with clumps or multiple nuclei (i.e., mergers) have high 
values, with M 20 > —1.1. In general, the further the 
brightest pixels are from the center of the source, the 
closer M 20 is to a value of zero. 

The right-hand panel of Figure shows the relation¬ 
ship between M 20 and the Gini coraficient for both the 
red and blue filters. For comparison we plot with crosses 
GANDELS galaxies selected to have 1.5 < z < 2 .5 and 
M*_ = IT^Mq and morphologically analyzed by Wang 
et al. (2012). To compare our sample to morphologically 
similar galaxies, we also plot with triangles a sample of 
73 local (z < 0.2 4) ULIRGS s t udied with HST in the rest 
frame o ptical bvlBorne et al. (2000). These sources were 
used by [Lotz et al. 


of z ~ 2 galaxies o 


To compare with their sample 


rserved with NICMOS, making them 
a suitable comparison set to our sources as well. Their 
morphological properties, including G, M 20 and A as well 
as a multiplic i ty cla ssification, are presented in Table 5 
of Lotz et al. (20041. The sources with “Double” nuclei 
are colored green. 

The dashed line indicates the separation be tween “nor- 
mal” and disrupted galaxies, as defined in [Lotz et 


(20041. All but one of our red quasars reside above the 
line m at least one of the filters. The local ULIRGs, 
particularly those with a “Double” nucleus morphology, 
have the lowest M 20 values and most closely approach 
the red quasar sample. 

However, the M 20 values of the red quasar sample are 
systematically shifted toward higher values with respect 
to the ULIRGs. Only the most extreme ULIRGs with 
“double” nuclei are consistent with the red quasar sam¬ 
ple. The higher M 20 values of the red quasar sample in¬ 
dicate host galaxy light distributions where the brightest 
pixels tend to be farther away from the central nucleus. 
We note that it is conventional to plot M 20 with lower 
values to the right, so the red quasars have the highest 
M 20 values and appear to the left in Figure [l5l 

The third metric that has been shown to effectively 
separate normal from merging galaxies is the rotational 
asymmetry, A. The standard prescription for comput¬ 
ing A involves rotating the image of a source, /, by 180° 
(creating the rotated image /iso) and producing asym- 
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Figure 15. Left - Gini coefficient versus E{B — V) measured in the red filter (top panel) and the blue filter (bottom panel). We see no 
correlation between the two quantities. Right - M 20 versus Gini coefficient for the red quasar host galaxies measured in the blue and red 
filters, and color coded accordingly. The two measurements for each object are connected by a dotted line. Plotte d for comparison with 
crosses are the same quantities computed for dusty star forming ga laxies at z 2 from CANDELS (|Wang et al.|2012|| and the circled crosses 
are X-ray sources in that sample. Local (z < 0.24) ULIRGs from|Borne et al.|||2000l are plotted with triangles and ULIRGs i dentified as 
having “Double” nuclei are colored green. The dashed line is the separation between “normal” and disrupted galaxies defined by |Lotz et al.| 
(|2004||. Our red quasar hosts are consistent with the same Gini values as the comparison samples, but only the most extreme ULIRGs with 
“'double” nuclei have M 20 values consistent with the red quasars. 


metry maps by taking the difference between /igo and 
I. The asymmetry maps are then used in the following 
formalism to quantify asymmetry: 


A = 


EiJ U(iJ) - -^18o(hj)l 


(9) 


Since our sources are dominated by a strong point source 
at the center, we experimented with rotating and sub¬ 
tracting the original image, as well as the PSF-subtracted 
frame (second column of Figure and found that the 
latter prod uces cleaner PSF-removed asymmetry maps. 
Figure [T6| shows the asymmetry images, I — Iiso for all 
eleven quasars as well as the PSF for the F160W images. 
In addition to clearly exposing the underlying disrupted 
host galaxies for most of the sources, we can re-examine 
three sources previously deemed undisturbed. F2M0943 
continues to show no underlying structure, F2M2222 ex¬ 
poses a small jet-like protrusion near its core that was 
hidden by the residual PSF light in the second and fourth 
panels of Figure F2M1359, which we characterize as 
being serendipitously reddened by the intervening galaxy, 
shows no asymmetry in its image, consistent with its host 
having a smooth symmetric profile as found by Galfit. 

We compute A for the 12 red quasars following equa¬ 
tion 1^ using the images shown in Figure To avoid 
PSF-related effects to affect our computation we exclude 
the innermost 16 x 16 pixels from our analysis and include 
only unmasked regions. The asymmetries we find have 
very high values {A = 0.99 — 1.6), s ignificantly highe r 
than those found for local ULIRGs in Lotz et al. (2004). 
However, we also found that the calculation of A is ex¬ 
tremely sensitive to how the background is defined as 
well as whether segmentation maps are used, versus the 
full image. 

Gollectively, the Gini coefficient, M 20 and A all suggest 


a high degree of asymmetry and dumpiness, as expected 
from merger-induced tidal effects and star formation. Ta¬ 
ble lists the Petrosian radii, Gini coefficients and M 20 
computed for the two filters. 

7. DISCUSSION OF INDIVIDUAL QUASARS 

In this section, we discuss each source individually in 
increasing redshift order, as they appear in Figuresj^and 

and note unique aspects of the morphological fits and 
host galaxy properties that we can infer. Throughout 
this section, statements about the “projected distance” 
to a companion system implicitly assume that the com¬ 
panion is at the redshift of the quasar. 

7.1. F2M1341 

The PSF-subtracted residual frame in Figure shows 
smooth, arc-like emission above and below the central 
point source. These are fit by two Sersic components 
in both filters. The southern component has a disk-like 
Sersic index of npiosw = 1-3 with an effective radius of 
Re = 6.4 kpc at a projected distance of 11.4 kpc. The 
redder component’s Sersic index and effective radius both 
have larger fitted values but are highly uncertain. The 
projected distance of the red component is a more distant 
12.4 kpc. 

More intriguing is the northern component which has 
well-separated blue and red emission concentrations. 
The blue component is nearer to the quasar, and has 
a de Vaucouleurs-like profile with npiosw = 4.8 while 
the red component is fit by a shallow, Gaussian or disk¬ 
like profile with npieow = 0.64. The two components 
are separated by ~ 4 kpc. In addition, there are two 
faint, red point-like sources to the east which we did not 
include in the fit as they may be faint low-mass Galactic 
stars. 
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Figure 16. Asymmetry maps for the F160W images produced by rotating the masked PSF-subtracted im^e, 7, by 180° to produce the 
rotated image, Iiso, and taking the difference (7 — 7i8o)- We use these maps in the numerator of equation!^ to compute A. The bottom 
right panel shows the asymmetry of the F160W PSF. These maps clearly show the complex and detailed struSure surrounding the residual 
of the point source in most of our sources; the exceptions are F2M1359 and F2M0943. 


7 . 2 . F2M1359 

This source has a host galaxy component detected in 
both F105W and F160W bands whose Sersic indices 
indicate a significant bulge component, although the 
precise value of n is uncertain (npiosw = 2.54 ± 0.4, 
?^Fi 60 W = 8.0 ± 4.9). Since this source is the only object 
with a unambiguous detection of a single host galaxy 
component at the quasar’s position with no additional 
components, the most likely scenario is one in which 
there is no merger. Additionally, there is a large, lower- 
redshift galaxy to the north of the quasar whose Sersic 
index indicates that it is disk dominated (n = 1.35), sug¬ 
gesting that the reddening in this case is due to extinction 
from the extended disk of the larger nearby galaxy. 

This would mean F2M1359 is an accidental red quasar. 
With an extinction corrected absolute-AT-band magni¬ 
tude of —29.7 (AB) it has the eighth highest luminos¬ 
ity of our sample, which is still remarkably luminous as 
compared with the unreddened quasars (black po ints in 
Figure [l I . Am ong the 13 red quasars studied by Urru- 
tia et ah (2008), one source (F2M0834-I-3506) was found 
to also be a normal quasar reddened by an interven¬ 
ing galaxy, so statistically it appears that < 10% of red 
quasars may be due to reddening that is not intrinsic to 
the quasars’ host galaxies. 


7 . 3 . F2M1036 

A highly asymmetric system with strikingly separated 
blue and red components that are best fit by two separate 
Sersic components in each band. All the components 


have small Sersic indices (n < 1). Some red light is seen 
in the residual image north-east of the source, which is 
not fully captured by the htting routine and is likely 
responsible for the higher value in the F160W image 
(6.00 versus 2.92). 


7 . 4 . F2M13U 

Despite being near a large intervening galaxy, unlike 
F2M1359 (( |7.2[ ) this source shows independent evidence 
for a merging host galaxy. The multi-component residual 
seen in the PSF-subtracted frame (second column) to the 
north of the PSF is, on its own, suggestive of a merger. 
There is also a faint point-like source in the model for 
the companion galaxy, which can be better seen in the 
inset of Figure 17 that shows a 6" x 6" image of this 
quasar in the F105W filter. This source was better ht 
by a Sersic profile with the large and unphysical indices 
■'T'Fiosw = 20.00 and n-Fieow = 7.17. Such a profile is 
centrally concentrated and indicates that this source may 
be a second AGN or perhaps a luminous, compact clump 
of star formation. 

The combined optical plus infrared spectrum of the 
quasar shows complicated absorption, i ncluding blue- 
shifted absorption in Mg 1 1 and Fe II (Urrutia et al. 
2009 Glikman et al. 2012), as seen in the optical-to- 
near-infrared spectrum shown in Figure [TtI While some 
of the UV absorption may be due to dust in the inter¬ 
vening galaxy, there are clearly in situ absorbing systems 
indicative of a merger. 
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Figure 17. Optical through near-infrared spectrum of 
F2M1344-h2839 demonstrating its unusual spectra l shape and 
classification as a FeLoBAL by|Urrutia et al.| (|2009||. The broad 
absorption features in Mg II and te 11 are marked with a red line 
and extend for 2000-5000 km s“^. The inset shows a 6" x 6" 
image in the F105W filter showing extended emission suggestive 
of tidal tails or a disrupted companion galaxy. There is also a 
faint point source visible in the frame. 

7.5. F2M0921 


This source has one of the clearest indications merging 
galaxies offset from the position of the quasar. There is 
also a point source 4" away that we ht in this model. 
From the Galfit photometry, the color of this second 
point source in the two filters is —0.11. When we cor¬ 
rect this color for the offsets to the UKIDSS passbands 
that were derived in Section the point source has 

Y — H = —0.33. The top panm of Figure [T8| plots the 

Y — H color of our source compared to colors of quasars 
as a function of redshift, based on the synt hetic UKIDSS 


colors derived for quasars (Tables 25 — 27 of Hewett et al. 


20061. The filled circle is the source’s Galfit colors. At 


z ~ 1.8 (the redshift of this source is z = 1.791) the 
Y — H color ranges between —0.1 and 0.2 mags. The 
point source source is significantly bluer than a typical 
unobscured quasar at this redshift. We therefore rule out 
this source as a companion quasar, which would be at a 
projected distance of 35 kpc away. 

The bottom panel of Figure 18 shows the Y — H color 


bles 13 and 14 

Hewett et al. 

(2006). The colors of this 

source are consistent with a w 

nite dwarf. And using the 


absolute magnitudes corresponding to the temperatures 
best-agreeing with this source’s color places such a white 
dwarf between ^ 40 and 330 pc. 


7.6. F2M0738 

This source is ^ 0.5"away from a companion galaxy, 
which - if at the same redshift as the quasar - is at a pro¬ 
jected distance of 16.4 kpc from the quasar. The galaxy 
is well fit by reasonable Sersic indices (npiosw = 2.50 
and npieow = 3.73) and corresponding effective radii of 
4.3 kpc and 3.2 kpc, respectively. We also fit an un¬ 
derlying host with Sersic indices of npiosw = 1-85 and 
■raFieow = 2.77, consistent with a bulge plus disk hybrid. 

7.7. F2M1427 

This system has one of the most complex morpholo¬ 
gies of our sample. The raw image shows a protruding 
structure to the south-east side of the quasar, and the 
PSF-subtracted frame in Figure 15 shows complex struc¬ 
ture, extended diffuse emission and ^ 5 red point sources. 



2000 4000 6000 8000 lOOOC 
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Figure 18. Predicted Y — H colors of quasars (top panel) and 
white dwarfs (botto m panel) from the UKIDSS colors derived by 
|Hewett et al.| (|2006|l are plotted to compare with the point source 
seen ~ 4" away from red quasar F2M0921 (black circle in the top 
panel; horizontal line in the bottom panel). The color of this source 
is too blue to be a quasar at this redshift, but has colors consistent 
with a Galactic white dwarf. 


The radio contours are circular and symmetrically cen¬ 
tered on the optical peak with no evidence of elongation 
lined up with the optical feature. However, the FIRST 
beam has a FWHM of 5"with 1"8 pixels, and higher res¬ 
olution radio images may reveal more complex structure. 

We model the significant emission by three Sersic com¬ 
ponents which capture much of the flux, but are likely 
not physically representative of the host systems. We 
conclude this because the F105W and F160W fit pa¬ 
rameters do not agree well with each other. This is 
most likely another example of a very comple x, multi- 
component merger similar to the systems seen in Urrutia| 


et al. (2008), but with a loss of the low surface brightness 
features needed to fully reconstruct the details of system. 


7.8. F2M0943 

This source would not converge with any physically 
meaningful added Sersic component. A single PSF com¬ 
ponent yielded a fit with reduced values of 51.7 and 
73.8 in F105W and F160W, respectively. When a Sersic 
component is added, the reduced improves greatly 
(though still leaving behind a strong residual) to 15.6 
and 23.1 in the two filters, respectively, but the effec¬ 
tive radii for the Sersic components are unphysical, at 
0.01 and 0.03 pixels, respective!}!^ In addition, Galfit 
assigns the added Sersic component more flux than the 
PSF component (i.e., the Sersic component is brighter 
than the PSF component by 0.7 — 0.8 magnitudes). 

We consider possible explanations for this poor fitting 
outcome. One possibility is that this source, because of 

Adding a second PSF rather than a Sersic component per¬ 
forms more poorly, yielding reduced values of 21.2 and 24.9 in 
the blue and red filters, respectively. 
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its higher redshift, was observed over two orbits. Breath¬ 
ing of the telescope or imperfect image combining via 
astrodrizzle may have affected the shape of the PSF 
in this field. We investigated the shapes of the three 
stars from the full reduced imaged of F2M0943 that went 
into creating our master PSF in both filters and did not 
find significant systematic offsets in the FWHM or other 
shape parameters as compared with stars in our other 
fields. We also constructed ‘mini-PSFs’ out of the stars 
in each field to look for evidence of differences between 
this field’s PSF compared with the others, and did not 
find any differences. The asymmetry map for this source 
also shows the strongest residuals in the point source, 
and does not show an obvious underlying host. 

Ruling out large PSF variations, we interpret the prob¬ 
lematic nature of fitting this source as either (a) there is 
something intrinsically different about this source (i.e., 
it is a luminous quasar hosted by an extremely luminous 
and centrally concentrated host), or (b) this object, being 
the most luminous quasar in our sample, is so bright that 
the Poisson noise dominates the residuals in the PSF fit¬ 
ting. Given that the residual image has low asymmetry 
and that the interpretation in (a) is highly unphysical, 
we propose that (b) is the more likely explanation. In 
support of this interpretation, we note that this source 
has the highest dynamic range (as defined in Section]^ 
and that residuals of this sign ificance are seen in some of 
the most luminous sources in Floyd et al. (2004). 


7.9. F2M2222 


This source is well-fit by a point source with the only 
nearby galaxy seen at a projected distance of 32 kpc. The 
nearby galaxy is a blue, disk-dominated (nFi 25 w = 1-62 
and npieow = 1.64) galaxy with an effective radius of 
Re = 6.1 kpc and Re = 5.8 kpc in the F125W and 
F160W filters, respectively, assuming it is at the same 
redshift. The residual image shows excess blue clumps in 
this galaxy suggestive of enhanced star formation. How¬ 
ever, it is not clear whether this galaxy is physically as¬ 
sociated with the quasar and its large projected distance 
indicates that this quasar is likely not hosted by an ob¬ 
vious merger. 


7.10. F2M1531 

This object is well-modeled by a central source plus 
a host galaxy that is 2.0 and 1.1 magnitudes fainter in 
F125W and F160W, respectively, and has a Sersic index 
indicative of a bulge (nFi 25 W = 3-64 and npieow = 4.48). 
However, Figurej^shows that some residual PSF flux re¬ 
mains in the image after the point source is subtracted. 
While there is a clear detection of an extended host 
galaxy, we consider its fitted flux to be an upper limit. 

Additionally, an elongated feature to the south is fit by 
a Gaussian profile (nFi 25 w = 0.67 and upigow = 0.61) 
at a projected distance of 12.3 kpc. This may be a tidal 
tail or the core of an interacting companion galaxy. 


7.11. UKFS0030 

This source is the lowest-luminosity and highest- 
redshift source in this sample (de-reddened absolute K- 
band magnitude of —30.99 at z = 2.299). It is the only 
red quasar o bserved from the UKFS sample of |Glikman| 


et al. (2013). This source is well-fit by a Sersic com¬ 
ponent atHie quasar location (with a central position 


at a projected distance of 0.9 kpc from the quasar) but 
only in the F125W filter. The component parameters 
are physically consistent with a concentrated bulge/disk 
hybrid (upiosw = 2.74) and Rg = 1.5 kpc, which may be 
interpreted as a nuclear starburst. There is also a nearby 
component to the northwest, at a projected distance of 
^ llkpc. 


8. MERGER STATISTICS FOR RED QUASARS 

Using a wide range of surveys from the literature for 
which AGN host morphologies are studied, Treister et al. 
(2012) examined the merger fraction among AGNs as a 


function of luminosity and redshift and found a strong de¬ 
pendence on luminosity over three orders of magnitude. 
The merger fractions ranged from 4% in Swift/HAT- 
detected AGN ( largely low redshift systems at z < 0.05, 
Koss et al. 2011) with Lboi ^ 10"^^'^ erg s~\ up to 85% for 
the red quasars studied by|Urrutia et al.| (|2008|), whose 
bolometric luminosities were the highest in the studied 
sample (Lboi 10^®'^ erg s“^). 

Here we place the z ~ 2 red quasars on th e merger 
fraction vs. l uminosity plot shown in Figure 1 of |Treistei^ 
et al.| (|2012D. Out of the eleven quasars studied in this 
paper, at inost three (F2M1359, F2M0943, F2M2222) 
sources do not show evidence for a galaxy merger. How¬ 
ever, we concluded that F2M1359 is not an intrinsically 
dust-reddened quasar and is likely obscured by an in¬ 
tervening galaxy lying along the line of sight, reducing 
the denominator in our calculation to ten intrinsically 
reddened quasars. F2M2222 is a more ambiguous case, 
leaving just one quasar (F2M0943) as showing no evi¬ 
dence at all for a merger. Therefore, we conservatively 
compute a merger fraction of 8/10 or 80%, although it 
could be as high as 90%. 

To compute the bolometric luminosities of the z ~ 
2 red quasars, we use bolome t ric co rrection from the 
quasar SED of Richards et al. (2006). Because of red¬ 
dening, we cannot use the optical or near-infrared mag¬ 
nitudes that we have present ed throughout this w ork. 
For the thirteen red quasars in Urrutia et al. (2012), full 
SED modeling was performed including Spitzer ItVi spec¬ 
tra and MIPS photometry out to 160 /xm, to compute 
their bolometric luminosities. Without data spanning 
such a broad wavelength range, we estimate the bolo¬ 
metric luminosities of our quasars by matching them to 
the Wide-F ield Infrared Space Explorer (WISE; [Wright 
et al. 2010) all-sky source catalog. All of our sources 
are detected within 1" in all four WISE bands. We use 
the longest-wavelength band, W4, whose effective wave¬ 
length is 22.0883 fim (corresponding to rest-frame wave¬ 
length between 6.1 /xm and 8.1 /xm), in order to mini¬ 
mize the effects of dust extinction and probe the intrinsic 
quasar emission. A t these wavelengths, t he bolometric 
corrections from the Richards et al. (2006) SEEp^ are all 
a factor of ^ 8. Using this method, our quasare lumi¬ 
nosities have a range log(Lboi) = 47.8 — 48.3 (erg s“^). 

Figure shows the merger fraction in AGN samples 
across many orders of magnitude in AGN luminosity. 
Gircles and triangles are data from a variety of AGN 
host galaxy studies and was incorporated into a meta- 
analysis of AGN triggering mechanisms by jTreister et al.j 

We use the SED that is made up of all the SDSS quasars in 
that sample. 
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(2012). This study of merging red quasar hosts at z ~ 2 
adds the most luminous AGN sample to this plot (red 
star). Our z ^ 2 red quasars are significantly more lu¬ 
minous but their merger fraction is commensurate with 
the 85% found at z = 0.7 (red circle). The dotted and 
dot-dot-dot-dashed lines are the parametrized li near and 


logarithmic fits to the data (equations 1 and 2 of Treister 
etahjlM^ . 


Tor a physical interpretation, w e compare ou r resul t 
to the variability-driven model of Hickox et al. (2014). 
Blue and red lines represent the p redictions tor merger 


fraction of AGN at z = 0.75 (the Urrutia et al. 
sample) and z = 2 (this work), respectively. 


2008 

While 


the error bars on our data point are too large to favor 
any particular model definitively, we note that while the 
lower-luminosity points agree better with the model that 
includes “mergers, interactions, and irregulars” (solid 
line), at the luminosities of our sources, the merger frac¬ 
tion is in better agreement with the model that only in¬ 
cludes “mergers and interactions” (dashed line). The 
latter model effectively excludes the role of minor merg¬ 
ers, which is a reasonable consideration at the luminosity 
regime of our sample. Therefore, it is likely that the role 
of minor mergers in black hole accretion declines with lu¬ 
minosity and is yet another lever that must be calibrated 
in simulations of cosmic AGN fueling. 

With a surface brightness limit of ~ 24 mag arcsec“^ 
(Figure^ the depth of our imaging is sensitive to major 
mergers out not sensitive to minor mergers. Thus the 
merger fraction we report in this paper is effectively that 
due to major mergers only. The fraction is consistent 
with the predicted contribution of m ajor mergers hosting 


the most luminous AGN according to Hickox et al. (2014 
dashed lines in our Figure [I^. The solid lines in Figure 
w hich include minor mergers as predicted by |Hick^ 
et al.| suggest that were our images deep enough to detect 
minor mergers, 100% of the sample would show evidence 
of minor or major merging activity. 


McLure et al. 

(19991), Kukula et al. 

(2001) 

Dunlop et al. 

(2003) and F 

Loyd et al. (2UU4) were also ot blue, lumi- 


their hosts are largely passive elliptical galaxies and they 
find no evidence for mergin g hosts. The most l uminous 
blue quasars presented by Floyd et al. (2004), reach¬ 
ing My < —28 at z < 0.5, also have elliptical profiles 
with only one (out of 17) objects showing clear signs of a 
merger. This is consis tent with the idea th at the quasars 
presented here and in Urrutia et al. (2008) were selected 
to be reddened Type I (i.e., showing broad emission lines) 
and are therefore in an earlier stage of the merger-driven 
evolutionary sequence, before the host galaxies’ mor¬ 
phologies settle into a virialized elliptical galaxy. 

This implies that even more heavily reddened quasars 
than the F2M and UKFS samples should show even 
stronger signs of mergers. Such sources may be hard 
to detect, but could be found by selecting a sample of X- 
ray sources with comparable redshifts and luminosities 
with no optical counterpart or spectroscopic evidence for 
AGN activity. Since such sources are expected to be rare, 
a large-area survey would be needed. S uch surveys are 
now bec oming available, e.g.. Stripe 82X dLaMassa et al. 


2013b|a) and the XMM -LSS surveys ( Fierre et al.||2l)L)7 


Ghiappetti et al.||2013 ). 



Figure 19. Merger fraction a s a function of bolom etric luminosity 
using the data from Table 1 of|Treister et aL|(|2012|| (black points). 
The red star represents the sources m tins work. 1 he dotted and 
dot-dot-dot-dashed black lines re present the linear an d logarith¬ 
mic fits to the points presented in |Treister et ah ]i |MT2l ) . Our data 
have effectively as high a merger fraction (?5Uyo) as tne next highest 
point (red circle; 85%) of [Urrutia et al]] | |2008^ but are more than 
an order of magnitude more luminous. 1 he blue and red lines show 
the p redicted merger fract ion based on the variability-based model 
from [Hickox et al.j (|2014|| realized at 2 : = 0.75 and z = 2. The 
solid line includes ‘‘mergers, interactions, and irregulars” in deter¬ 
mining the fraction, while the dashed line only considers “mergers 
and interactions”, which effectively ignores minor mergers. The 
merger fraction in our sample agrees well with the “mergers and 
interactions” model at 2 = 2 and suggests that major mergers play 
the most significant role in fueling the most luminous quasars at 
2 = 2 . 


Interestingly, the morphologies of our quasars suggest 
that they are in an early stage of the mergers, with clear 
independent companions rather than asymmetr ic or dis¬ 


rupted single hosts. H owever, the models (e.g., Hopkins 
et al.|2005 2006 2008) predict that the AGN are bright¬ 


est tor longest in the late merger stage. One possible 
explanation is that these red quasars are associated with 
the early merger stages when the dust is typically far 
from the nucleus, rather than the late-stage merger when 
the dust can be more concentrated around the nucleus 


(as seen in t he heavily obscured Type H quasars in Lacy 
|et al. I [20071 which seem to be mostly late-stage merg- 
ers). Another interpretation is that the red quasars are 
tracing group interactions in which the quasar host has 
already been through a merger, but additional mergers 
are ongoing. 


9. SUMMARY AND CONCLUSION 

We have conducted a morphological study with the 
WFC2/IR camera on HST of eleven z ^ 2 dust reddened 
quasars that are intrinsically the most luminous systems 
at this redshift. This is the highest redshift sample of 
HST-imaged dust-reddened quasars. Visual inspection 
shows clear evidence of mergers in at least eight of the 
eleven systems; only one source does not appear to reside 
in a merging host. 

We performed careful PSF modeling in the three broad 
near-infrared filters in order to remove the quasar light 
and study the morphologies of their underlying host 
galaxies. We modeled the galaxies with parametric fit- 
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ting to a Sersic profile and see significant merging com¬ 
ponents in eight systems. One system proved to be red¬ 
dened by an intervening galaxy and is not technically 
a “red quasar”. We use the Sersic profiles to integrate 
the total flux of the host galaxies to study their proper¬ 
ties and find that these galaxies are more luminous than 
the massive galaxies at z 2 studied by the CANDELS 
survey. Although redder than the overall systems’ col¬ 
ors, the red quasar host galaxies are mostly bluer than 
the CANDELS galaxies. Three sources are significantly 
redder, suggesting a very dusty star bursting host. We 
interpret the diversity of colors coupled with high lumi¬ 
nosities as consistent with these systems still exhibiting 
ULIRG-like properties as a result of a major merger. 

We conduct a non-parametric analysis of the PSE- 
subtracted images, measuring their Gini parameter, M 20 
and asymmetry and find broad consistency with local 
ULIRG properties, particularly with the most extreme 
double-nuclei ULIRGs, further supporting them being 
hosted by major mergers. 

Compared with studies of the merger fraction in AGN 
at different luminosities, our sample is more luminous 
by two orders of magnitude than the previous most- 
luminous sample of red quasars studied in this manner 


[mil 

(jUrrutia et aL|[2008|) and has a comparable merger frac¬ 
tion! When added to AGN samples from other stud¬ 
ies at a variety of AGN luminosities on a plot of host 


galaxy mer ger fraction versus AGN luminosity (Treister 
et al. 2012|) and co mpared with variab ility-driven black 
hole growth models ( Hickox et al.||20i4 ), our data favors 
a model in which black hole growth in the most lumi¬ 
nous AGN at z ^ 2 is dominated by major mergers. 
We conclude from this study that major mergers are the 
dominant drivers of black hole growth at the highest lu¬ 
minosities as far back as the epoch of peak quasar and 
star formation activity. 


TOP CAT (Taylor 2005) and an OS X widget form 


of the JavaScript C osmology Calculator (Wright 2006 


Simpson et al.|[2013l were used while preparing this pa¬ 


per. 
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Table 1 

Table 1: High Redshift Red Quasars Imaged with HST WFC3/IR 


Name 

(1) 

R.A. 

(J2000) 

(2) 

Deck 

(J2000) 

(3) 

J 

(mag) 

(4) 

Ks 

(mag) 

(5) 

(mag) 

(6) 

Redshift 

(7) 

Orb 

(8) 

Filter 

(9) 

Exp 

(sec) 

(10) 

(mag/s/”2) 

(11) 

Filter 

(12) 

Exp 

(sec) 

(13) 

(mag s“^ 
(14) 

F2M1341-I-3301 

13:41:08.11 

+33:01:10.2 

17.83 

16.81 

1.00 

1.715 

1 

F105W 

897 

23.8 

F160W 

1597 

24.0 

F2M1359-I-3157 

13:59:41.18 

+31:57:40.5 

17.86 

16.72 

0.88 

1.724 

1 

F105W 

897 

23.6 

F160W 

1597 

23.9 

F2M1036-I-2828 

10:36:33.54 

+28:28:21.6 

17.93 

17.15 

0.85 

1.762 

1 

F105W 

897 

23.6 

F160W 

1597 

24.0 

F2M1344-I-2839 

13:44:08.31 

+28:39:32.0 

17.48 

16.66 

0.30 

1.770 

1 

F105W 

897 

24.0 

F160W 

1597 

24.0 

F2M0921-I-1918 

09:21:45.69 

+19:18:12.6 

17.70 

16.48 

1.20 

1.791 

1 

F105W 

897 

23.3 

F160W 

1597 

23.4 

F2M0738-I-2750 

07:38:20.10 

+27:50:45.5 

17.99 

17.18 

0.99 

1.985 

1 

F105W 

897 

23.9 

F160W 

1597 

24.0 

F2M1427-I-3723 

14:27:44.34 

+37:23:37.5 

18.09 

16.99 

0.60 

2.168 

1 

F105W 

997 

23.9 

F160W 

1597 

23.9 

F2M0943-I-5417 

09:43:17.68 

+54:17:05.5 

16.97 

16.15 

0.30 

2.232 

2 

F105W 

2194 

24.2 

F160W 

3193 

24.2 

F2M2222-0202 

22:22:52.79 

-02:02:57.4 

18.39 

17.07 

0.77 

2.252 

2 

F125W 

1794 

24.2 

F160M 

3194 

24.3 

F2M1531-I-2423 

15:31:50.47 

+24:23:17.6 

17.58 

16.60 

0.75 

2.287 

2 

F125W 

1794 

24.2 

F160W 

3194 

24.3 

UKFS0030-I-0025 

00:30:04.96 

+00:25:01.3 

19.32 

18.01 

0.90 

2.299 

2 

F125W 

1794 

23.9 

F160W 

2794 

24.0 


^ The 3 it surface brightness limit quoted here is a per pixel value (0.06” X 0.06”) 


Table 2 

Point Spread Function Parameters 


Filter 

Number of Fields 
Archive Program 

-Nstars 

FWHM 

(") 

F105W 

7 

8 

35 

0.147 

F125W 

9 

3 

29 

0.2094 

F160W 

9 

11 

46 

0.1782 
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Table 3 

Galfit Multi-Component Fitting Parameters 


Source 

Name 

Model 

Type 

Y 

(mag) 

Re^ 

(pixels) 

n 

Dist 

(") 

X?ed 

H 

(mag) 

Re^ 

(pixels) 

n 

Dist 

(") 

X?ed 

F2M1341 

PSF'^ 

18.59±0.25 



0 

4.777 

17.3±0.01 



0 

5.430 

Comp 

Sersic 

23.24±0.13 

12.7±2.3 

1.31±0.26 

1.35 


22.3±1.2 

61.1±144.8 

7.64±6.10 

1.44 


Comp 

Sersic 

21.69±0.57 

79.7±69.5 

4.8±1.6 

1.59 


22.4±0.4 

12.6±0.7 

0.64±0.08 

2.06 


F2M1359 

PSF*> 

20.74±0.31 



0 

2.407 

18.48±0.31 



0 

13.505 

Host 

Sersic 

21.51±0.31 

3.5±5.8 

2.54±0.40 

0.11 


19.6±1.4 

425±1015 

8.0±4.9 

0.19 


F2M1036 

PSF'’ 

18.58±0.01 



0 

2.917 

17.97±0.01 



0 

6.003 

Comp 

Sersic 

22.23±0.01 

5.55±0.06 

0.30±0.03 

1.46 


21.24±0.01 

8.20±0.06 

0.20±0.02 

1.38 


Comp 

Sersic 

21.84±0.01 

7.22±0.07 

0.13±0.01 

0.81 


22.75±0.06 

11.0±0.7 

0.92±0.14 

0.77 


F2M1344 

PSFb 

18.42±0.26 



0 

5.652 

18.39±0.16 



0 

11.061 

Host 

Sersic 

22.14±0.39 

9.2±6.0 

1.03±0.71 

0.54 


19.25±0.55 

4±220 

0.01'' 

0.02 


Comp 

Sersic 

19.6±1.6 

3202±12685 

20.00'' 

1.21 


20.13±0.67 

233±269 

7.17±2.25 

1.25 


Comp 

Sersic 

22.85±0.04 

6.64±0.16 

O.lOiO.05 

1.58 


22.15±0.02 

7.7±5.1 

0.04'' 

1.52 


F2M0921 

PSF'’ 

19.13±0.16 



0 

3.905 

17.54±0.16 



0 

9.269 

Host 

Sersic 

20.16±0.31 

4±62 

o.oG 

0.06 


19.47±0.32 

6.1±2.3 

19.3'' 

0.31 


Comp 

Sersic 

21.98±0.07 

lO.lil.l 

2.43±0.26 

1.23 


19.79±0.67 

193±271 

10.9'' 

1.45 


Comp 

Sersic 

18.69±1.84 

3570±11295 

11.74'' 

2.10 


20.1±1.5 

408±1325 

12.3'' 

2.33 


AGN? 

PSF 

21.27±0.01 



4.13 


21.38±0.18 



4.11 


Comp 

Sersic 






21.12±0.29 

23±26 

11.7'' 

4.15 


F2M0738 

PSF'’ 

18.63±0.26 



0 

3.287 

17.50±0.27 



0 

3.047 

Host 

Sersic 

22.22±0.39 

16.9±6.3 

1.85±0.73 

0.55 


22.92±0.80 

3±6 

2.8±1.4 

0.46 


Comp 

Sersic 

21.36±0.02 

7.02±0.17 

2.13±0.06 

1.86 


20.65±0.01 

6.60±0.14 

3.93±0.11 

1.84 


F2M1427 

PSF 

18.56±0.01 



0 

6.084 

17.70±0.01 



0 

5.220 

Comp 

Sersic 

21.82±0.08 

15.25±2.17 

3.68±0.38 

2.26 


19.95±0.43 

287±274 

11.86'' 

2.27 


Comp 

Sersic 

24.36±0.08 

1.91±0.34 

1.79±0.84 

1.38 


23.61±0.11 

2.1±0.8 

6±3 

1.39 


Comp 

Sersic 

23.11±0.12 

9.7±1.7 

1.78±0.33 

1.54 


21.22±0.23 

40±15 

4.6±0.8 

1.60 


F2M0943 

PSF” 

17.57±0.31 



0 

15.564 

16.85±0.01 



0 

23.063 

F2M2222 

PSF 

18.40±0.01 



0 

2.889 

17.97±0.01 



0 

9.699 

Comp 

Sersic 

21.56±0.02 

12.32±0.34 

1.62±0.05 

3.93 


21.53±0.03 

11.65±0.52 

1.64±0.08 

3.93 


F2M1531 

PSF'’ 

18.06±0.16 



0 

3.912 

17.75±0.30 



0 

7.646 

Host 

Sersic 

20.01±0.15 

4.4±5.2 

3.64±0.39 

0.19 


18.86±0.31 

1.9±5.9 

4.48±0.33 

0.12 


Comp 

Sersic 

24.17±0.04 

5.09±0.33 

0.6±0.2 

1.55 


23.93±0.04 

5.4±0.4 

0.61±0.25 

1.54 


UKFS0030 

PSF 

19.55±0.31 



0 

1.997 

18.72±0.01 



0 

4.238 

Host 

Sersic 

21.11±0.30 

2.9±5.9 

2.74±0.36 

0.11 







Comp 

Sersic 

22.52±0.01 

5.48±0.11 

0.65±0.05 

1.35 


22.38±0.01 

5.83±0.13 

0.71±0.06 

1.30 



^ The distance to the Sersic radius in pixels, which can be converted to arcseconds using the image scale of 0.060”/pixel. 


The PSF magnitude reprted here is comprised of the integrated flux from two PSF components, converted to a magnitude 
using equation^ 

The PSF magnitude reprted here is comprised of the integrated flux from PSF component plus a Sersic component that is 
sharply concentrated at the position of the PSF, with unphysical parameters. The total PSF magnitude is computed with 
equation 

These parameters are flagged by Galfit as being outside the range of acceptable values, however, the fit resulted in an 
acceptable Xred enabling a capture of the total flux in the host components. We do not report errors for these parameters. 
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Table 4 

Quasar Magnitudes and Colors 

SExtractor MAGJIUTO GALFIT primary PSF'^ GALFIT combined PSF*’ 


Name 

U 

(mag) 

V 

(mag) 

U-V 
(mag) 

U 

(mag) 

V 

(mag) 

U-V 

(mag) 

U 

(mag) 

V 

(mag) 

U-V 

(mag) 

E{B - V) 
(mag) 

F2M1341 

18.66±0.07 

17.38±0.02 

1.28±0.02 

19.24±0.25 

17.29±0.01 

1.95±0.04 

18.59±0.25 

17.29±0.01 

1.30±0.06 

0.57 

F2M1359 

20.51±0.06 

18.27±0.03 

2.24±0.03 

20.99±0.31 

18.50±0.31 

2.49±0.28 

20.74±0.31 

18.48±0.31 

2.26±0.28 

0.50 

F2M1036 

18.67±0.03 

18.11±0.02 

0.55±0.02 

19.12±0.01 

18.51±0.01 

o.eiio.oi 

18.58±0.01 

17.97±0.01 

o.eiio.oi 

0.47 

F2M1344 

18.41±0.03 

17.78±0.02 

0.63±0.02 

19.02±0.26 

18.50±0.16 

0.52±0.14 

18.42±0.26 

18.39±0.16 

0.03±0.15 

0.07 

F2M0921 

19.01±0.03 

17.62±0.02 

1.40±0.02 

19.34±0.16 

17.60±0.16 

1.74±0.14 

19.13±0.16 

17.54±0.16 

1.59±0.13 

0.65 

F2M0738 

18.64±0.02 

17.63±0.02 

1.01±0.02 

19.25±0.26 

17.68±0.27 

1.57±0.22 

18.63±0.26 

17.50±0.27 

1.13±0.21 

0.49 

F2M1427 

18.58±0.02 

17.85±0.02 

0.73±0.01 

18.56±0.01 

17.70±0.01 

0.86±0.01 

18.56±0.01 

17.70±0.01 

0.86±0.01 

0.27 

F2M0943 

17.70±0.02 

16.91±0.07 

0.78±0.05 

18.75±0.31 

18.10±0.01 

0.65±0.11 

17.57±0.31 

16.85±0.01 

0.72±0.10 

0.11 

F2M2222 

18.51±0.02 

17.99±0.02 

0.52±0.01 

18.40±0.01 

17.97±0.01 

0.43±0.01 

18.40±0.01 

17.97±0.01 

0.43±0.01 

0.33 

F2M1531 

18.00±0.02 

17.50±0.02 

O.SliO.Ol 

18.68±0.16 

17.96±0.30 

0.72±0.20 

18.06±0.16 

17.75±0.30 

0.30±0.18 

0.32 

UKFS0030 

19.71±0.09 

18.87±0.03 

0.84±0.04 

19.55±0.31 

18.72±0.01 

0.83±0.10 

19.55±0.31 

18.72±0.01 

0.83±0.10 

0.32 


Note. — The magnitudes presented in this table are directly measured from the HST blue (F105W or F125W) and red (F160W) 
bandpasses as pseudo rest-frame U and V magnitudes, respectively. 

^ In sources where two PSFs were used to fit the quasar flux, we consider the brighter fitted component to be the primary component. 
These magnitudes include integrated quasar magnitudes arising from fitting a single PSF, two PSFs, or a PSF plus Sersic components. 


Table 5 

Non-Parametric Measures for Red Quasar Host Galaxies 
Object Blue Filter Red Filter 


Name 

-^Pet 

(arcsec) 

Gini 

Af20 

-^Pet 

(arcsec) 

Gini 

M 20 

E(B - V) 
(mag) 

F2M1341 

1.55 

0.61 

-0.39 

0.44 

0.43 

-0.39 

0.57 

F2M1359 

2.12 

0.48 

-0.60 

2.12 

0.49 

-0.61 

0.50 

F2M1036 

0.41 

0.22 

-0.73 

0.43 

0.25 

-0.78 

0.47 

F2M1344 

0.36 

0.37 

-0.47 

0.51 

0.57 

-0.55 

0.07 

F2M0921 

1.76 

0.61 

-0.38 

2.05 

0.62 

-0.35 

0.65 

F2M0738 

2.11 

0.64 

-0.66 

0.43 

0.31 

-0.76 

0.49 

F2M1427 

0.43 

0.50 

-0.47 

0.42 

0.20 

-0.60 

0.27 

F2M0943 

0.42 

0.20 

-0.52 

0.51 

0.59 

-0.27 

0.11 

F2M2222 

0.41 

0.52 

-0.33 

0.53 

0.54 

-0.26 

0.33 

F2M1531 

1.16 

0.71 

-0.26 

0.92 

0.60 

-0.33 

0.32 

UKFS0030 

1.66 

0.66 

-0.40 

0.42 

0.39 

-0.84 

0.32 



